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We find an expression for mass of spherically symmetric system solving Euler-Lagrangian equa¬ 
tion by Homotopy Perturbation Method. With the help of this expression and the Einstien Field 
equations we obtain an interior solutions set. Thereafter we explain different aspects of the solution 
describing the system in connection to mass, density, pressures, energy, stability, mass-radius ratio, 
compactness factor and surface redshift. This analysis shows that all the physical properties, in 
connection to brown dwarf stars, are valid with the observed features except that of stability of the 
model which seems suffers from instability. 
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I. INTRODUCTION 

There has always been search for the interior solution 
of spherically symmetric system. Considering different 
methods and approximations researchers tried to find out 
more than several hundreds different type of interior so¬ 
lutions. But very few solutions made its physical accep¬ 
tance in all aspects describing the system so far. 

However, in the present paper we have studied a spheri¬ 
cally symmetric system of radiating star under the homo¬ 
topy perturbation method (HPM). This is a series expan¬ 
sion method used in the solution of nonlinear partial dif¬ 
ferential equations, in the present case the Einstein field 
equations of general relativity. The method in principle 
employs a homotopy transform to generate a convergent 
series solution of differential equations. The HPM was 
introduced and developed by He [IH3 and others (sl-fT^. 

In one of his earlier works He Q proposed a cou¬ 
pling method of a homotopy technique and a perturba¬ 
tion technique to solve non-linear problems. The author 
argued that in contrast to the traditional perturbation 
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methods, the proposed method does not require a small 
parameter in the equation and the results reveal that the 
new method is very effective and simple. Later on several 
other researchers have employed the HPM in the diverse 
field of pure and applied mathematics (as an exhaustive 
appliance), and in physics and astrophysics (as a new 
field of application) to solve related non-linear differen¬ 
tial equations in an extraordinary simplified way [l3l - [l^ . 

In connection to self-gravitating radiation model 
Sorkin et al. [T^ have examined the entropy of self- 
gravitating radiation confined to a spherical box with 
finite radius in the context of general relativity. Their 
results are expected to supports the validity for self- 
gravitating systems of the Bekenstein upper limit on the 
entropy to energy ratio of material bodies. 

Rahaman et al. [T^ proposed and analyzed a model 
for the existence of strange stars where observed masses 
and radii are used to derive an interpolation formula 
for the mass as a function of the radial coordinate. 
The interpolation technique has been used to estimate 
the cubic polynomial that yield the following expres¬ 
sion for the mass as a function of the radial coordinate 
m(r) = ar^ — br^ + cr — d with a, b, c and d all being 
numerical constants. Their analysis is based on the MIT 
bag model and yields physically valid energy density, ra¬ 
dial and transverse pressures. 

However, in the above mentioned work of Rahaman 
et al. [l^ the target object was a strange star. In our 
present investigation we start with the intention to de- 
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velop a basic interior solution of the Einstein equations 
valid for any radiating model under a similar expression 
for the mass as a function of the radial coordinate m(r). 
Then we match our theoretically obtained solutions set 
with the observational results for practical validity of the 
model and find that our model is best fit for the brown 
dwarf star of EO type. 

Therefore, the motivation of our present work is to 
find the interior solution of spherically symmetric sys¬ 
tem considering radiation effects in a more general way. 
For this we consider a metric with unknown time-time 
component. By using the maximum entropy principle 
(MEP) we obtain an Euler-Lagrangian equation which 
gives second-order differential equation for mass of the 
spherical distribution. Then we find out solution for 
mass by the homotopy perturbation method (HPM). The 
Einstien field equations are solved for density, time-time 
component of metric and pressures of the system. Thus, 
basically taking data set arbitrarily to see the physical 
nature and from there to predict whether this type of ra¬ 
diation model can effectively be used is the main purpose 
of the present investigation. 

The work plan of the study is as follows: in Sec. II we 
develop the methodology for mass of the spherical sym¬ 
metric system of radiating star under (i) the maximum 
entropy principle, and (ii) the homotopy perturbation 
method. Next in the Sec. Ill we provide the Einstein 
field equations for self-gravitating radiation system and 
discuss several physical aspects of the solution describ¬ 
ing the system in connection to mass, density, pressures, 
energy, stability, mass-radius ratio, compactness factor 
and surface redshift. In Sec. IV it has been shown in de¬ 
tails that all the physical features in connection to brown 
dwarf stars of type EO are valid with the observed fea¬ 
tures except that of stability of the model which seems 
partially unstable. We have made some specific remarks 
on the radiation model under consideration in the Sec. 
V. 


II. THE METHODOLOGY FOR MASS OF THE 
SPHERICAL SYMMETRIC SYSTEM OF 
RADIATING STAR 


where, Pr and p are the radial pressure and density of the 
matter distribution. Compatible with spherically sym¬ 
metry, we assume the general energy-momentum tensor 
as 


= ip + - \p€ + ipt - ^p)v^Vu, (3) 

with = 1. 

Here = 1 and u^pi, = 0. The vector 

is the fluid 4-velocity of the local rest frame of the 
radiation and is the unit spacelike vector which is 

orthogonal to where the radial pressure pr is in the 

direction of 

For locally measured temperature T of black body ra¬ 
diation, the rest frame energy density and entropy den¬ 
sity respectively assume the forms as [l^ 

P = hT\ (4) 


s = (5) 

where & is a constant of order unity (in Planck units G = 
c = /i = fc=l), on the assumption that the number of 
species of radiation is of order unity . 

By substituting Eq. (4) into Eq. (5), one can easily 
obtain the following relation 

s = (6) 


where a = 

For the matter distribution up to radius r < R, the 
total entropy is given by 


rR 

S = Air s{r) 
Jo 

which can be simplified as 


1 - 


2m(r) 


1 -1 


r^dr, 


S = (47r)ia / Ldr, 
Jo 


( 7 ) 


( 8 ) 


A. The Maximum Entropy Principle 

We consider the metric of the spherical system 

ds^ = —gttir)dt^+^ - dr"^ +r^{dO^+sin^9d(p^), 

( 1 ) 

where m{r) is the mass distribution of the spherically 
symmetric body and gtt(r) is the time-time component 
of metric. These are functions of radial coordinate r only. 

The equation of state for the system with thermal ra¬ 
diation is given by 


where the Lagrangian L is given by 

( 9 ) 

with 

m'{r) = Airr'^p, (10) 

which can be shown from the Einstein field equations 

Gab=8nTab [ 3 . 

Also, we have Euler-Lagrangian equation in the form 


L={m'p 


1 - 


2m(r) 


Pr = ip, 


( 2 ) 


d 

dr 



dL 

dm' 


( 11 ) 









which reduces to 
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TO + m m — —mr — -= 0- (12) 

This is the master equation for a self-gravitating radia¬ 
tion system. If it is possible to solve this highly nonlinear 
differential equation, then one would get complete struc¬ 
ture of the self-gravitating radiation system and therefore 
could be used to describe a model of fluid sphere or stellar 
configuration. In the following section, we shall solve this 
equation by using a newly developed technique known as 
homotopy perturbation method (HPM) Q. 


B. The Homotopy Perturbation Method 

Primarily our task is to find the solution of Eq. dm) by 
using homotopy perturbation method. In this method we 
consider linear part of the differential equation as m' = 
dm/dr and the rest of the terms as non-linear. So the 
homotopy equation can be written as Q 


! / 7 / n ^ II 2 2 4toto' 

m—m^-\-ti mQ + mm “ -= b, 

(13) 

where h is a parameter and toq is the initial guess of the 
solution. For h = 0, we get initial approximation solution 
and for h = 1, one obtains the desired solution for mass. 

We assume the initial solution as toq = ar^. According 
to HPM, therefore, we consider the general solution as 

TO = Too -I- hmi -I- h^m2 -I- h^m^ + .... (14) 

Substituting this in Eq. (HU and equating coefficients 
of the different orders h, we get the solution for mass up 
to second order correction as follows 


36 


312 


m{r) = oTO -I- -^a TO -|- i^a'^r + SCior + C 2 , (15) 
0 oO 

where Ci and C 2 are constants of integration. 

The boundary conditions in the present system are 

to(0) = 0, m'{R) = 0, (16) 

with 

m'{r) = 3ar^ + 36a^r^ + -|- 6Ciar. (17) 
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So, we immediately obtain the constants of integration. 
Cl and C 2 in the following form 


Cl =- 


and 


R ^ o 52 2 „ 

— -b 6aR^ + —a^R 


C 2 = 0, 


(18) 


(19) 


where R is the radius of the spherical distribution. 
Therefore, substituting Ci in Eq. (HSI, we obtain 

/ \ 3 36 2 5 312 3 7 

m{r) = ar^ + ^ 


-3ar^ 


R 


52 


-b 6 ai?^ -b — 


III. THE EINSTEIN FIELD EQUATIONS FOR 
SELF-GRAVITATING RADIATION SYSTEM 

For the matter distribution given in Eq. (IlD, the Ein¬ 
stein field equations for the metric © can be written 
as 


2m' 

— = OTTp, 


2to _ ^1 _ ‘//m\ ^1 


Stt 


r J 9ttr 


( 21 ) 


( 22 ) 


- 1 - 


2 to 
r 

/to to 

V r 


2 5 tt 4 \gtt) 2.r gtt 


i + 14 

r 2 gtt 


+ (23) 


From Eq. m we get the density as 

m' (r) 


P = 


A'lrr'^ 


Therefore, substituting Eq. (HU into Eq. 
tain an expression for density as follows 


(24) 
we ob- 


P = 


47r 


2 2 312 3 4 

3a -b 36a^r^ H-o^TO 
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Again by using Eq. m and Eq. dSll), we solve Eq. 
for gtt as follows: 


gtt = K 


exp / 3(r-\m) 
r(r — 2 to)i/3 ’ 


(26) 


where A is a constant. 

It can be noted that for negative values of a and any 
value of r the exponential term becomes unity and hence 
we obtain 


K 


gtt 


r(r — 2 to)^/^ 


(27) 


Using boundary conditions for gtt, one can find out 
K, which ultimately gives the time-time component of 
metric as 


9tt = \ — 


^^ 4/3 _ 2rn^y/3 


r J (1 _ 2m(r) ^i/3 ’ 

where m{R) is a constant and can be given by 

m{R) = --aR^ - —a^R/' - —a^R\ 
2 5 7 


(28) 


( 20 ) 


(29) 
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IV. PHYSICAL FEATURES OF THE 
RADIATION MODEL 

A. Density and Mass 

To start with the basic physical studies preferably we 
choose parameter a = —1.73 x 10“^® and R — 0.065 so¬ 
lar radius and m{R) = 0.080 solar mass of EO class of 
brown dwarf stars (see the Link [l^). The reason to 
consider these data set lies on the physical background 
that our target is to investigate a radiating model which 
is compatible with brown dwarfs which are cool star-like 
objects that have insufficient mass to maintain stable nu¬ 
clear fusion in their core regions [IJ. It has been argued 
that although brown dwarfs are not stars, they are 
expected to form in the same way, and they should radi¬ 
ate a large fraction of their gravitational energy at near- 
infrared wavelengths. 

For the above specification we plot variation of mass 
and density as shown in Fig. 1 and 2. 




FIG. 2: Variation of density as a function of radial distance 
r for a = —1.73 x and R = 0.065 solar radius 



FIG. 3: Variation of radial pressure as a function of radial 
distance r for a = —1.73 x 10~^® and R — 0.065 solar radius 


FIG. 1: Variation of mass as a function of radial distance r 
for a = —1.73 x 10“^® and R = 0.065 solar radius 


B. Pressure and anisotropy 


From Figs. 3 and 4 we observe that the radial and 
tangential pressures and hence anisotropy (Fig. 5) are 
decreasing function of radial distance r. 


Using Eq. m and (I25p , we get an expression for radial 
pressure as 


Pr = 


1 

1^ 


3a -I- 36a^r^ -|- 

5 


--ff+6ai?3 + fa2i?5 

r V 2 5 


(30) 


whereas from Eq. (l23l) we get tangential pressure as 
follows 


1 


Pt = 


727rr^ 


„ / 2m 

3m r — 10m H- 

r 

(m' - si f 

41-rjL + 4 


(31) 


C. Stability 

In connection to stability analysis we have checked out 
the stability criteria given by Herrera [2lj. In this ap¬ 
proach (usually known as the concept of cracking or over¬ 
turning) the region for which the radial and tangential 
sound speeds < 0 and its numerical value less 

than unity is a stable region. We observe from the plot 
in Fig. 6 that for there is no change of sign but 

up to 41800 km its numerical value greater than unity. 
This implies that our configuration with EO type brown 
dwarf star of 0.065 Solar radius (i.e. 45227 km) is unfor¬ 
tunately suffers from instability within the region of core 
radius 41800 km. 

However, it is interesting to observe from the Fig. 7 
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FIG. 4: Variation of tangential pressure as a function of radial 
distance r for a = —1.73 x 10“^® and R — 0.065 solar radius 



FIG. 5: Anisotropy, A = {pt — Pr), of the pressures as a 
function of radial distance r for a = —1.73 x 10~^® and R — 
0.065 solar radius 


that outside the core our model is potentially stable one 
where clearly Vst - Vsr < 0 and its numerical value less 
than unity. Therefore, from the plot of Fig. 7 it is ob¬ 
vious that within the shell starting approximately from 
41800 km to 45227 km the value of sound speed differ¬ 
ences is less than one which is the condition of stability. 
Like wise the plot is maintaining same sign so there is no 
crack which is also another condition of stability. Hence 
only 7.57% of the shell (3427 km) near the surface the 
Brown Dwarf has a stable configuration. 


D. Energy conditions 

Let us now check whether all the energy conditions re¬ 
garding present radiating model are satisfied or not. For 
this purpose, we are considering the following inequalities 


FIG. 6: Sound speed ~ '^sr < 0 through out the region of 
core radius (41800 km) implying instability of configuration 
with EO type brown dwarf star of 45227 km where unit of r 
is in km 



FIG. 7: Sound speed v^t — v^r < 0 through out the region of 
the thin shell (3427 km) implying stability of configuration 
with EO type brown dwarf star of 45227 km where unit of r 
is in km 


and plotting graphs for each case: 

(i) NEC : p + Pr > 0, P + Pt >0, 


(ii) WEC '■ p + Pr > 0, p>0, p + Pt>0, 


{Hi) SEC '■ p + Pr > 0, P + Pr + ‘^Pt > 0. 

We see from Figs. 8 and 9 related to different energy 
conditions indicate that in our model all the energy con¬ 
ditions are satisfied through out the interior region. 
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E. Compactness and Redshift 


By definition compactness of a star can be given by 
u = m[r)lr which in our case takes the following form: 


1 - NEC - WEC^\ 

n 36 2 4 312 o fj „ / 

u = ar^-\ -- a^r^—Zar ( 

-h 6ai?^ -b 



5 35 ' 

. 2 

5 ; 


FIG. 8: The energy conditions applicable in the interior of the 
spherical distribution have been plotted against r for NEC 
and WECr 


(32) 



Fig. 10 shows that < |. 


Again, surface redshift is defined by the relation 


z = (1 - 2u)-^ - 1, (33) 


- WEC, -SEC 


FIG. 9: The energy conditions applicable in the interior of the 
spherical distribution have been plotted against r for SEC and 

WECt 


which in the present study turns out to be 


J 


z = - 1. (34) 

^^[1 - 2{ar2 + f a2r4 + - 3ar (f + 6aR^ + f 


The variation of surface redshift z is shown in Fig. 11 
which is physically within an acceptable profile. 


V. CONCLUSION 

In the present paper we have searched for an expres¬ 
sion of mass with spherically symmetric system to solve 
Euler-Lagrangian equation by Homotopy Perturbation 
Method. By employing this expression for mass and the 


Einstien Field equations we have obtained an interior so¬ 
lution. We have also explained different physical prop¬ 
erties of the solution describing the system. However, in 
the present model for radiating brown dwarf stars all the 
features are seen to be physically viable except that of 
stability of the matter distribution. 

In connection to these aspects with brown dwarf stars 
we have done a preliminary calculation taking the data 
set of Bhar et al. [2^ related to highly compact stars 
such as neutron stars and strange stars. This data set 
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FIG. 10: Variation of compactness as a function of radial 
coordinate a = —1.73 x 10“^® and R — 0.065 solar radius 


provides negative radial pressure and also in this case 
the model suffers from instability. 

Therefore, we suspect that - 

(1) radiation model may not be compatible with highly 
compact stars rather it is compatible with radiating 
brown dwarf stars. Actually, to relate radiation model 
with highly compact stars is also not justified due to 
their exhausted fuel system and hence burning as well 
as radiation stage which seems have been stopped much 
earlier. 

(2) instability may be inherent property of any radi¬ 
ating compact stars under HPM and MEP as observed 
in the present model with brown dwarf stars. Instability 
has also been observed in the extension of the work of 
Bhar et al. [2^ related to highly compact stars such as 
neutron stars and strange stars. 



FIG. 11: Variation of redshift 2 as a function of radial distance 
r for a = —1.73 x 10“^® and R = 0.065 solar radius 


As a final remark we would like to state that some 
more investigations are needed to perform with differ¬ 
ent methodology before coming to a definite decision re¬ 
garding (i) applicability of HPM and MEP to radiating 
and highly compact stars, and (ii) instability of radiating 
brown dwarf stars under HPM and MEP. 
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